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Abstract 
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I. INTRODUCTION 



Bose-Einstein condensates in a periodic optical lattice potential have attracted much in- 
terest. An optical lattice is an ideal system for studying a variety of solid state physics 
phenomena, such as Bloch oscillations [l|, Q|, Landau-Zener tunneling ^, Q|, and Josephson 
oscillations An optical lattice is an ideal crystal without lattice imperfections or impuri- 
ties, and parameters can be changed easily over a wide range. The experimental achievement 
of the superfluid-Mott insulator transition in an optical lattice p also proved the usefulness 
of optical lattices to understand basic properties of strongly correlated systems. Moreover, 
Bose condensates in an optical lattice are predicted to show novel phenomena such as dy- 



namical instability 



12, 



13, 



□ QDD 



111 ] and swallow-tail energy loops in the band structure 
141 ] . Recent studies have concentrated on the case of pure Bose condensates at zero 
temperature. The effect of a thermal cloud on the dynamics of the condensate excitations 
in an optical lattice have not been studied very much. This is the subject of the present 
paper. 

Recently, several experimental papers have reported results on the collective modes of 



Bose condensates in a one- dimensional optical lattice [15|, IrJ [l?! 13] and its damping at 
finite temperature [16]. In these experiments, a shift of the oscillation frequency in the 
presence of the optical lattice and a sharp change of the damping rate with increasing depth 
of the optical lattice have been observed. The measured frequency shift ^j] of the collective 
modes is in good agreement with the renormalized mass theory of Kramer et al. Q] . The 
damping of condensate oscillations in a ID optical lattice at T = has been measured jl7|. 
The damping of the condensate excitations at finite T has not been studied in any detail. 

Since our major interest in this paper is the thermal cloud of excitations, it is important 
to make a clear distinction at the outset between (a) The Bloch-Bogoliubov excitations asso- 
ciated with linearized fluctuations of an equilibrium Bose-condensate and (b) The stationary 
states of the time-independent Gross-Pitaevskii equation for the Bose order parameter. In 
a continuum model, the latter states can be described by the eigenfunction (we use a ID 
lattice for illustration) 

= e ikx u k (x), (1) 

where the condensate Bloch function satisfies the usual periodicity condition u^{x) = u k (x + 
Id), where d is the optical lattice spacing and I is an integer. Physically, <&1(x) corresponds 



to a solution of the GP equation with a superfluid flow in the periodic potential with the 
condensate quasi-momentum k. Recent theoretical literature (see, for example, Refs. 7] and 

has reported extensive studies of such condensate Bloch states, including their energy 
band structure and stability. The latter question can be studied by considering the dynamic 
fluctuations 8$k(x,q) around the equilibrium state with the generalized Bogoliubov 

excitation energy E q (k). These excitations are also described by a quasi-momentum q in 
the first Brillouin zone (BZ) and will be referred to as the Bloch-Bogoliubov excitations of 
the optical lattice. The thermal cloud of non-condensate atoms which is present at finite 
temperature is described as a gas of these Bloch-Bogoliubov excitations. As emphasized 
in the literature j^j], one must not confuse the energy bands of these Bloch-Bogoliubov 
excitations E q (k) with the condensate energy bands or Bloch eigenstates described by Eq. 

That is, we must distinguish between the "condensate" energy band and the "excitation" 
energy band. In our tight-binding model, the analogue of Eq. (0) is 

= e Md ^Hk), (2) 

where d is the optical lattice spacing and I is an integer. While one could generalize our 
analysis, we only consider the Bose condensate in the k = Bloch state, in which case Eq. 
(J2J) reduces to 

$° fc= o(0 = v 7 ^. (3) 

Here n c0 denotes the number of condensate atoms trapped in each well of the optical lattice, 
labelled by I. 

In this paper, we study the dynamics of Bose atoms in an optical lattice at finite temper- 
atures. We discuss a Bose-Hubbard tight-binding model using the Gross-Pitaevskii approx- 
imation, but generalized to include the presence of non-condensate atoms. We calculate the 
temperature dependence of the condensate atom number n c0 (T) in each lattice well using 
the static Popov approximation. This is needed in our calculation of the Landau damping of 
condensate modes due to coupling to thermal excitations. We also extend our earlier results 
for a !D optica! lattice Q to 2D aod 3D optica! lattices. A s „ in Re, Q, tor damp- 
ing processes to occur, the dispersion relation of the Bloch-Bogoliubov excitations E q must 
initially bend upward as the quasi-momentum q increases. This is referred to as "anoma- 
lous dispersion" and is also the source of 3-phonon damping of long wavelength phonons in 



superffuid 4 He |22j, [23]. This condition leads to a dramatic disappearance of all damping 
processes of phonon modes in a D-dimensional optical lattice when a = Un c0 /J > 6D, 
where U is the on-site interaction and J is the hopping matrix element. 

This paper is organized as follows. In Sec. Owe introduce the well-known Bose-Hubbard 
tight-binding model which describes Bose gases in an optical lattice. In Sec. Illl[ we dis- 
cuss the generalized discrete Gross-Pitaevskii equation and the characteristic changes in the 
Bloch-Bogoliubov excitation dispersion relation as a function of the dimensionless interac- 
tion parameter a = Un c0 /J. In Sec. IIV1 we introduce the static Popov approximation for 
optical lattices, and calculate the condensate fraction and a as a function of the optical 
lattice depth and temperature. In Sec. El we calculate the Landau damping of Bogoliubov 
excitations in an optical lattice of dimension D. We also briefly remark on Beliaev damping 
into two excitations. We made some concluding remarks in Sec. IVII 

As mentioned, some of results of this paper in a ID optical lattice were briefly described 
in Ref. j^lj. In the present paper, we give a more detailed derivation and discussion, as well 
as new results for 2D and 3D optical lattices. 

II. MODEL 

We consider bosonic atoms in an optical lattice potential 

D 

V op {r) = sE R J2^ 2 (k Xl ), (4) 

1=1 

where s is the usual dimensionless parameter describing the optical lattice depth in units 
of the photon recoil energy Er = h 2 k 2 /2m. D is the dimension of the optical lattice and 
d = ? = | is the lattice period. We only consider simple cubic lattices considered in recent 
experiments (], I^J. We call attention to the recent technique 1?J of producing a two- 
dimensional array of long, tightly confined condensate tubes by loading a Bose condensate 
into a deep 2D optical lattice potential, which prevents atoms from hopping between different 
tubes. With an additional ID optical lattice potential along a tube, an ideal ID system 
can be experimentally realized [l7|. One can also have an ideal 2D system by loading a 
condensate into a deep ID optical lattice and a shallow 2D optical lattice. We assume this 
experimental setup for the realization of ID and 2D optical lattices in the present paper. 
We also assume that the laser intensity determining the depth of the optical lattice wells is 



large enough to make the atomic wave functions well localized on the individual sites (i.e., 
where we can use the tight-binding approximation). The energy gap between the first and 
the second excitation bands is large compared to the thermal energy (2/c B T /E R <C s), and 
thus only the first band is thermally occupied. 

Within a tight-bindin g appro ximation, the Hamiltonian is effectively described by the 



Bose-Hubbard model [25, 



26, 27 



p. 

2' 



as 



H = - J ( a ) a i + a \ a j) + (^) 

<3,l> 3 

where a,j and a] are destruction and creation operators of atoms on the j-th lattice site, 
(j, I) represents nearest neighbor pairs of lattice sites. The first term describes the kinetic 
energy due to the hopping of atoms between sites. The hopping matrix element J is given 
by 

h 2 V 2 



J 



drw*(r) 



2m 



+ V C 



opV 



(6) 



where Wj(r) is a wave function localized on the j-th lattice site, and m is the atomic mass. 
Expanding the optical lattice potential around the minima of the potential wells, the well 
trap frequency is lu s = s 1 ^ 2 ^-. Approximating the localized function as the ground state 
wave function of a harmonic oscillator with frequency u s at the potential minima of j-th 
site 



one obtains 



J 



,% 2 s 



( 



mw,\ D / 4 



E 



R 



V 7rh 



,1/2 



exp 



--2F (r - r '' 



;S (l + exp (_ s -l/2)) 



(7) 



(8) 



Here s 1 ^ 2 ~ y/ 2m (^ E R) ^ can ^ g interpreted as a WKB factor for tunneling in an optical lattice 
potential which has height sEr and well width d. 

The second term in Eq. J5J) describes the interaction between atoms when they are at 
the same site. We assume that atoms can move along z direction in ID case and in xy plane 
in 2D case. The on-site interaction U depends on the dimensionality of the optical lattice. 
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U is given by 26] 



U = g drj0 ± (r ± )| 4 / d*K-(*)| 4 , 



(ID) 



9 j d ^|0||(^)| J dxdy\wj(x,y)\ , (2D) 
g J dv\ Wj (v)\\ (3D) 



(9) 
(10) 

(11) 



where g = Anh 2 a/m and a is the s-wave scattering length. Here 0j_(r_|_) = 

(^nt) 1 ^ 2 ex P (~^ft r2 ±) an< l ^11 ( z ) = (~^h~) 1 ^ ex P i~~lh~ z2 ) are ^ ne g roun d state wave func- 
tions in optical lattice well traps for confining atoms in 1 and 2 dimensions. Approximating 
the localized function simple Gaussian, one obtains 



U_ 
E~r 



9 



2 3 / 2 ad 1 



(2n) 3 / 2 a 2 ± a s ^l 2 a\ 
g 2 3 / 2 a 



(27r) 3 / 2 a||a2 tt 1 / 2 ^ 
g 2 3 /V/2 a 



(2nfl 2 al 



d 





(ID) 


s 1/2 , 


(2D) 


s 3 '\ 


(3D) 



(12) 
(13) 
(14) 



where a s = * I , a± 



n 



muij_ ' 



and an = \ I — — . 



III. THE GENERALIZED GROSS-PITAEVSKII EQUATION 

In order to investigate the collective excitations of a Bose condensate in an optical lattice, 
we restrict ourselves to the superfluid phase of the system described by Eq. ((SJ). We start 
from the Heisenberg equation of motion for a,j(t) (we set h — 1 from now on), 



i_ a .(*) = [ aj {t),H] 



-J^2ai(t) + Ua](t)aj(t)aj(t), 



(15) 
(16) 



where (I) means that I runs over the nearest neighbor sites of site j. In the presence of Bose 
condensation, we write aj(t) = &j{t) +ipj(t), where $?(£) = (cij(t)) is the condensate wave 
function at site j and ipj{t) is the non-condensate field operator (with (ipj{t)) = 0). This 
non-condensate operator ipj satisfies the usual Bose commutation relations, 



o. 



(17) 
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The equation for the condensate wave function is obtained by taking an average of the 
Heisenberg equation. If, as usual, we neglect the anomalous correlations rrij(t) = (ipjipj) and 
(ipfyj^j) , we are left with a "generalized" form of the discrete Gross-Pitaevskii equation 



which includes the mean field of non-condensate atoms 



2< 



i^.( t ) = -jJ^Mt) + U( n p) + 2^(t))$,(t). (18) 

(i) 

Here, n^(i) is the number of condensate atoms on the j-th lattice site and hj(t) = (V>!'0j) is 
the number of the non-condensate atoms on j-th lattice site. The time-dependent Hartree- 
Fock mean field 2Uhj{t) in Eq. ()18j) arises from the thermal gas of non-condensate atoms. 
Equation (fTH|) reduces to the usual discrete Gross-Pitaevskii equation [^, 0, |3| when all the 
atoms are assumed to be in the condensate (i.e., if we set hj = 0). 

Introducing the usual phase and amplitude variables for the condensate, %(t) = 
n c( i ) e w J (t) j Eq is equivalent to 

dn c At) 



dt 



2J E \A#kw Mm - Oj(t)), (19) 



d0;(*) _ T ^ l n i( t \ C os{6 l {t)-6j{t))-U{n c j {t)+2hj{t)) 




dt ^ V r^(t) 

(20) 

Equation (fT9"j) is the continuity equation for the condensate and Eq. ()20|) often called the 
Josephson equation. The energy of a condensate atom £j(t) reduces to the equilibrium 
chemical potential fi C Q of the condensate in static thermal equilibrium (8j, n c - and hj are 
independent of site index j) 

fi c0 = -zJ + U(n c0 (T) + 2h°(T)). (21) 

Here z is the number of the nearest neighbor sites, n c0 is the number of the condensate 
atoms per site in equilibrium, and h° is the number of the non-condensate atoms per site 
in equilibrium. The solution of Eq. (j2U|) in static equilibrium is 6o{t) = — A*cO^- From 
Eq. (|19J). one finds that the Josephson current between the j-th and Z-th lattice sites is 



J 3 {t) = 2J^n%t)n c l {t) sin(0,(*) - 9j(t)). 



The well-known Bogoliubov excitation spectrum for a uniform optical lattice is easily 
obtained from the above equations 22,0, 3]- Considering small fluctuations from equi- 
librium, n c j(t) = n c0 + Sn^t), 6j{t) = o (t) + 56j(t) and hj(t) = h° + <%(*), Eqs. (HHJ) and 
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FIG. 1: The Bogoliubov excitation spectrum E„ in a ID optical lattice plotted as a function of the 
quasi-momentum q in the first Brillouin zone. The spectrum is plotted for a < 6 and a > 6. The 
dotted lines give the low q phonon dispersion relation E q = cq. 



PH) reduce to 



d5n%t) 
dt 

dt 



-2Jn c0 ^(^(t) -66j(t)) 



J 



(Snm - Sn%t)) - U (Sn%t) + Mntf)) . 



(22) 
(23) 



Ignoring the non-condensate atom term (5fij = 0), the solution of these coupled equations 
are the normal modes 



50j{t) = 56(^y^^- E ^\ 5rf 3 (t) = 5n c (q)e i[q - r ^ E ^. 
The Bloch-Bogoliubov excitation energy in an optical lattice is given by 

£ q = ^/e°(e°+2tfn<*), 

where 

i=i 

For small q, this spectrum is phonon-like E„ ~ eg, with the phonon velocity 



(24) 



(25) 



(26) 



c = V2Jd 2 Un c0 
8 



Un c0 



rrr 



(27) 



where = is an effective mass of atoms in the first excitation band of the optical 



lattice [19. l3l| . This T = excitation spectrum in ID case is shown in Fig. 1, for two values 
of the dimensionless interaction parameter a = Un c0 /J. 

We call attention to an important feature of the dispersion relation E q in Fig. ^ con- 
sidered as a function of the ratio a. For a < 6, the excitation energy E q bends up before 
bending over, as q approaches the BZ boundary. This behavior is analo gou s to the so called 



22 



23j,|32j. For a > 6, 



"anomalous dispersion" of the phonon spectrum in superfluid 4 He 
in contrast, the spectrum simply bends over as one leaves the low q (phonon) region. This 
feature will play a crucial role when we consider damping processes in optical lattices. The 
excitation spectrum in 2D and 3D optical lattices also exhibit this kind of spectrum. How- 
ever, the critical value of a then depends on the direction of q, since simple cubic optical 
lattices in 2D and 3D do not have rotational symmetry. The crucial effect of this anomalous 
dispersion on damping processes will be discussed in detail in section |V| 



IV. STATIC POPOV APPROXIMATION 

In this section, we calculate the condensate fraction n c0 / n and the parameter a = U n c0 / J 
as a function of the optical lattice depth s and the temperature. This is needed as input 
to the calculations in Sec. We base our discussion on the static Popov approximation 

Substituting aj = $j + ipj to K = H — fiN, where N = a ] a j is the total number of 



9 



atoms, one obtains, 



K = K + K 1 + K 2 + K 3 + Ki, (28) 

k = -jE( $ ^ + $ ^)-^£i $ / + y£i $ /> ( 29 ) 

(hi) 3 3 



(30) 



3 V (0 J 

K * = - J £ + - ^ £ fyk 

m 3 

+^ E {^]f + + ($*) 2 ^ 2 ) , (31) 

j 

K 3 = UJ2 + > ( 32 ) 

i 

Here, we use the Hartree-Fock-Bogoliubov-Popov approximation that takes into account the 
third and fourth order terms of t/>] within a mean- field approximation, but neglects terms 
involving the anomalous averages (ipjipj) and (w0/0j), namely we use 

K 3 ^ 2*7 ^ n,- ($*^- + , K,~2UY, njtyfa ■ (34) 

i i 

Within this approximation, we have 

= K + K[ + K' 2 (35) 

+ u ( nC j + 2 ^;) i>3 (36) 
v (i) J 

K> i = £($*+#&) -**£$& 



3 

U 



+ 2 E + 4 ^M- + (*i) 2 #) > ( 37 ) 

3 

where rij = n^ + n^ = |$j| 2 + (ipjipj). The coefficients of the linear terms in ■0, in Eq. (J36j) 
identically vanish using the fact that satisfies the generalized Gross-Pitaevskii equation 
Eq. IjlBjl . As a result, = 0. 



10 



As discussed in Sec. HI we consider a Bose condensate in static thermal equilibrium in 
the k = Bloch state, and hence $j = \/n c0 . Using the value of the condensate chemical 
potential in thermal equilibrium given in Eq. (j21j) . the remaining term K' 2 in Eq. (J3HJ) 

reduces to 



£(n^j) 2 +n c0 $). (38) 



07} 

17 
2~ 





We next introduce Fourier components, 

^ = 4wH a ^' r ^ ( 39 ) 

where the number of lattice sites in one direction is denoted by I and hence the total number 
of sites is I D . The Momentum sum is over the first Brilloiuin zone of the optical lattice. 
Substituting Eq. ((HSJ) into Eq.flHBJ), one finds 

K' 2 = E( e 2 + Un^a^ + — EK a -q + a Q a -q)- ( 4 °) 

q^O q^O 

We can diagonalize K = Kq + K' 2 by the Bogoliubov transformation j^ . |3| , 

Oq Wqflq -|- ?7 q Ci_ q , 

«-q = V - q fl q + U-qOtq- ( 41 ) 

If we assume that a q and a q obey the usual Bose commutation relations, we obtain the 
following conditions for w q and i> q , 

| Mq | 2 - | Wq | 2 = 1. (42) 

From the condition for the diagonalization of K, we find that w q and u q have to satisfy the 
following equation, 

Un c \u\ + v\) - 2(e q + Urf^u^ = 0. (43) 

Solving both Eq. (J42|) and Eq. (|43|) . one can easily derive the parameters for the Bogoliubov 
transformation (we take u q and t> q real), 

«H(t4"-KH- (44) 
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We have here introduced the Hartree-Fock (HF) excitation spectrum 



£ q = [e (q) + 2U(n c0 + n )}-ii c0 
,2/5< , rr„cO 



4Jsin 2 (y) +Un c ". (45) 



Putting all these results together, our total Hamiltonian has been reduced to 



q^O q^O 



= -^V° + 4n°) + \ ]T (E q - £ q ) + £ £ q a q a q . (46) 

q^O q^O 

Here iV c o = n c0 I D is the total number of condensate atoms. The Bogoliubov-Popov excita- 
tion spectrum E^ appearing in Eq. (j4^j) is identical to Eq. (j2HJ), except that now n c0 (T) 
is the temperature-dependent number of condensate atoms in any given lattice well. The 
Bogoliubov-Popov excitation spectrum can be obtained by ignoring the non-condensate fluc- 
tuation (5hj = 0) but using the temperature-dependent number of condensate atoms n c0 (T) 
as given by Eq. (j22j) and Eq. ([23)1 . We thus assume that the thermal cloud is always 
in static equilibrium when dealing with the time-dependent density fluctuations Sn^it) of 



cloud has been successfully used 
which we refer the reader for 



the condensate. The approximation of a static therma 
in other discussions of condensate collective modes j3, 
further discussion. 

It is easy to verify that £" q in Eq. (|25|) reduces to the HF energy in the limit J ^> U n c0 , 
i.e., when aCl. This HF limit corresponds to setting the Bogoliubov amplitudes w q = 1, 
i>q = in Eq. (jUj) . In dealing with Bose gases trapped in harmonic potentials, one can 
always use this HF approximation j^] for the excitations describing the thermal cloud as 
long as the kinetic energy of the atoms (~ k-e,T) is much larger than the interaction energy 
(Un c0 ). In contrast, apart from the limiting case of a <C 1, we must always use the full 
Bogoliubov spectrum E^ to describe the thermal cloud composed of excitations in the first 
band of an optical lattice. 



Expressing n°(T) in terms of these Bogoliubov-Popov excitations, we have (27 

1 

1 



28, 



q^O 



q^O 

12 



where f°(E q ) = [exp(f3E q ) — is the usual Bose distribution function. The number of 
condensate atoms rf° at a site is found by solving Eq. (}4"Tj) self-consistently for a fixed value 
of the total site density n. The condensate fraction n c0 In in a .D- dimensional optical lattice 
is shown in Fig. |21 We take n = 2 and use the parameters from Ref. [17] The spurious finite 
jump in the condensate atom number n c0 at the transition temperature T c is an inherent 
problem of the Bogoliubov theory in a uniform gas (for further discussion see, for example, 
Ref. I3J). 

Strictly speaking, there are no solutions of Eq. ()47|) at finite temperature for an infinite 
optical lattice in ID and 2D because of the divergent contribution from excitations with 
small momentum, in accordance with the well-known Mermin-Wagner-Hohenberg theorem 
36^ . However, for the finite systems discussed in this paper, this divergence is not present. 
That is, Eq. ()47|) has a solution describing a finite value of the condensate n c0 (T) below T c . 

In Fig. El we plot the parameter a = Un c0 / J as a function of the temperature for several 
values of the optical depth s. The dimensionless interaction parameter a and the results in 
Fig. El will play an important role in the discussion in the rest of this paper. Since we limit 
our discussion to the first energy band of the optical lattice, our results only apply when 
s ^> 2kftT/En. Higher excitation bands would be thermally populated if we consider lower 
values of s and would have to be considered. 



V. DAMPING OF BOGOLIUBOV EXCITATIONS 

In this section, we discuss the damping of condensate modes in optical lattices. In Sec. 
IV Al we derive a formal expression for the damping of Bogoliubov excitations using the 
thermal Green's function formalism. In Sec. IV Bl we use this expression to calculate the 
Landau damping of Bogoliubov excitations in ID lattice. Section IV CI extends this analysis 
to 2D and 3D optical lattices. 



A. Green's function technique 



In this section, we take = K + K% + K 2 in Eq. (|33|) as the bare zeroth-order 



Hamiltonian and treat K* perturbation (We neg 
use the simplest Bogoliubov theory to diagonalize K-q 



ect the fourth order term K4). We 
28l | . We neglect the fourth order 



27 
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FIG. 2: The condensate fraction n c0 /n in D-dimensional optical lattice as a function of tempera- 
ture. In this paper, the height of the optical lattice potential (in units of £r) is always denoted by 
s. 14 
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FIG. 3: The demensionless interaction parameter a = Un c0 (T)/ J, plotted as a function of temper- 
ature, for several values of the optical well depth s. The number of condensate atoms at a lattice 
site n c0 (T) is given in Fig. |21 The dashed line shows the critical value a c = 6D, above which there 
is no damping. 
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term in the fluctuation K±. The chemical potential for Kb is given by /ibo — —zJ + Un c0 
(from the condition K\ = 0) . K-i can be diagonalized by the Bogoliubov transformation Eq. 
(jUJ) and one obtains Eq. ()46p. 

Using Eq. (J39j) and Eq. (|41j) . K 3 can be re- written in terms of the Bogoliubov quasipar- 
ticles 



-^3 2 ^ ] J ^pi,P2;P3 ( a p 3 a pi a P2 ~l~ ^pa^pi^r^) 

pi,P2,P3^0 

+ ^ ] -^pi,p2;P3 ^ a p3 Q; -pi Cl; -p2 a -P3 a Pi a P2^ 5 (^) 



Pl,P2,P37^0 

where the matrix elements are 

/ jyCO 

M piip2;p3 = ZUjj^Y^ [i U Pl U P3 + V P 3 - W P1 M P3) M P2 

G 

(49) 

/ ^c0 

^pi,p 2 ;P3 = U\l —jj ("^Pl^^PS ~~ n pi n P2 W P3) ^pi+P2,P3+G) (50) 

G 

and G is a reciprocal lattice vector. The Kronecker delta <5p 1+P2i p 3+ G expresses the conser- 
vation of quasi-momentum of the three excitation scattering processes (Umklapp processes 
are associated with G^O). 

The thermal Green's function for Bogoliubov quasiparticles is defined by (see e.g. 

G q (n - r 2 ) = (T T {a^a^n))), (51) 

where 

a q (r) ee e^ T a q e-^ r , a q (r) ee e KT a\e~ Kr . (52) 

The angular brackets (...) indicates taking an average for thermal equilibrium and T T gives 
the usual time ordered product of operators. Fourier transformation of this imaginary time 
Green's function is given by 

G^T) = ±J2 G ^ e ~ iUJnT > ( 53 ) 

where ui n is the boson Matsubara frequency u n ee ^ (n is an integer). The zero-th order 
Green's function is 

G l(^n) = : , 1 F , (54) 
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P2>^n - 




q, u n 



FIG. 4: The second order self-energys. 



where _E q is the Bloch-Bogoliubov excitation energy given in Eq. (1251) . 

The second-order self-energy terms are given by the usual diagrams shown in Fig. 4. 
Fig. |U (a) describes a condensate excitation of (quasi) momentum q being absorbed by 
an excitation p 2 of the optical lattice thermal gas, lending to a thermal excitation with 
momentum pi. This describes Landau damping processes Fig. |U (b) describes 

a condensate excitation of momentum q decay into two excitations with momentums pi and 



p 2 . This describes Beliaev damping 
of the Bogoliubov excitation 



4l|. These diagrams give the second-order self-energy 



E q (iu; n ) 



£^n) 



lUi - lUJ ri 



E 2 

Pl,P2^0 LU l 



(55) 



The summation over Matsubara frequencies is calculated by the well-known formula (see 
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Ref. |22|) 

^2gM = ~fdzg(z)f(z), (56) 

where f°(z) is the Bose distribution function and the contour encircles the poles of g(z). 
Using Eq. (|5S|). one obtains 

The damping of excitations is given by the imaginary part of the self-energy 

r q = -Im£ q (£ q + i5), (59) 
where 5 — > +0. This gives the damping of Bogoliubov excitations in an optical lattice, 

r q = rj + r;, (eo) 

^ = vr J2 I^P^H/V^pJ-f^pJl^^q-^+^pJ, (6i) 

Pl,P2#0 

^ = \ E iMp^Hl + r^pJ+Z^pJlM^-^-^pJ- (62) 

Landau damping is expected to be dominant at higher temperatures where there are 
thermally excited quasiparticles. In contrast, Beliaev damping Tq is due to a decay process 
and can arise in the absence of thermally-excited excitations. Beliaev damping is possible 
even at T = [i.e. f°(E p ) = 0] and is expected to be dominant at low temperatures. 



B. Landau damping in a ID optical lattice 

In this subsection, we calculate the Landau damping given by Eq. (|61|) for a ID optical 
lattice. For this, the energy conservation condition E n + E p = £q +p+ G = -Eq+ P needs to 
be satisfied, where G is a reciprocal lattice vector. The solution of the energy conservation 
condition E q + E p = E q+P for a ID optical lattice is illustrated in Fig. 03 42]. First, we draw 
an excitation spectrum and mark q on the line. Then we draw the spectrum line again with q 
as origin. If those two lines intersect, solutions consistent with energy conservation condition 
for three excitation process exist. The intersection can be taken as (q + p, E q + E p = E q+P ). 
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q q + p 

FIG. 5: The Bogoliubov excitation energy E q in a ID optical lattice, for a < 6. The intersection 
of the two dispersion curves at p + q shows that the energy conservation condition is satisfied in 
this case. 



Clearly, the condition for these two dispersion curves to have intersections requires that the 
dispersion relation E q first bends up as q increases, before bending over. From Fig. 1, we 
see that the ID optical lattice dispersion relation E q has this feature only for a < 6 and thus 
Landau damping can occur only when a < 6. If the intersection (q + p,E q + E p ) is outside 
of the first Brillouin zone, it has to be reduced in the first Brillouin zone by subtracting a 
reciprocal lattice vector G n = ^p- {n is an integer), corresponding to an Umklapp process 

42|- 

The values of (q,p) satisfying the energy conservation condition for ID optical lattice are 
shown in Fig. |3 For a given value of q, we see that as a — > 6, the value of p decreases 
to zero. There is no solution for a > 6, indicating the disappearance of the damping for 
an excitation E q with any value of q. In Fig. 5, the curves of the solution of the energy 
conservation condition never cross the dashed line at q + p = ir/d. Thus Umklapp scattering 
processes (G n ^ 0) do not contribute to Landau damping in our present single-band model. 

When the excitation of wavevector q has a wavelength much larger than the thermal 
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FIG. 6: The values (q,p) satisfying the energy conservation condition E q + E v = E q+p , for several 
values of the interaction parameter a. 

excitation p (i.e., when q <^ p, it/d), one finds that = c from the energy conservation 
condition E q + E p = E q+p in Eq. (JoTj). The Landau damping of the excitation E q = cq 
comes from absorbing a thermal excitation E p with a group velocity equal to c j^]. This 
condition requires 

dE P E v A1A ■ P d P d 
— — = — 4 J a sin — cos — 

dp E p 2 2 

= c = V2Jd 2 Un c0 . (63) 

This equation can be written in the form 

16X 2 + 8(a - 2)X + (a 2 - 6a) = 0, (64) 

where X = sin 2 ^. The unique value of p sucn that in Eq. (J53j) equals c is found to 
be given by the condition 



• 2 / PoA ~(a-2) + y/2^ + 2) 

sm ( -y J = 4 • ( 65 ) 
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This expression is only valid for a smaller than 6, such that po ^> q. 
For q <C p, where we can use the following approximations, 



! Un c0 



u q — 



2cp V 2Un c0 J ' 



Jsin(pd) ([/n c0 ) 2 , , . 

_z Mp — 9^, (68) 

Jsin(prf) (f/n c0 ) 2 , 

<W *v p - 2 ^ jgg g rf » ( 69 ) 

the matrix element in Eq. (jfil)) reduces to 

Mq^q+p = 2U d — — [ {UpUq+p + VpV q +p — V p U q +p) U q — (u p Uq+ p + V p V q + p — U p V q +p) V ? J 

orr /"""V 6 ? C' 100 (Un<<>) 2 Jd . A r^f~ ,_. 

From the energy conservation condition 

d£p = £p 2Jdsin d (n) 

dp £ P v ; 

the matrix element can be approximated by 




E P \( q V /2 Ko 



1/2 



M «,^lF^K^j (72) 



where N c0 = n c0 I. This is the analogue of Eq. (10.76) in Ref. |43j. Using the delta function 
for energy conservation to integrate over p, the Landau damping is finally reduced to 



\ 2 

,l a f E poY I e °o , E po \ PJUqd 



q 8(a + 2) \e° J \E Po + E m ) sinh 2 ^ ' 



2 



One finds there is no damping when a > 6 (see Fig. 3 for a as a function of s and T). In 
fact, diverges at a = 6, but this is due to the ID nature of our system. This divergence 
does not occur in 2D and 3D lattices. 



C. Landau damping in 2D and 3D optical lattices 

We next discuss the energy conservation condition in a 2D optical lattice. We imagine 
a surface in a three dimensional space which satisfies the Bogoliubov dispersion relation 
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(q x ,q y , Eq). Then, we draw a new dispersion Eq in the three dimensional space, with a 
point on the surface (qi x , qi y , E m ) as origin. That is, we draw (q x , q y , -E q - qi + E m ). If those 
two surfaces intersect, the energy conservation condition Eq 1 + Eq 2 = E m+Ci2 is satisfied, the 
intersection being given by (qi x + <l2x, Qix + Q2y, E ni+n2 ). For this condition to be satisfied, the 
surface (q x , q y , Eq) has to be above the other surface (q x , q y , Eq_q 1 +Eq) around (qi x , qi y , Eq^. 
Since the Bogoliubov spectrum is phonon like Eq ~ cq for small q, the maximum gradient 
of Eq at (qi x ,qi y , EqJ must be greater than c. Since |V q -Eq| q=qi is the maximum gradient 
of Eq at qi, this condition is equivalent to the requirement: 

E I 

I V q .E q | q = qi = 2Jd-^-dsm 2 q lx d + sin 2 q ly d > c. (74) 

Equation (fT^j) is the 2D version of the condition for the Bogoliubov spectrum of a ID 
optical lattice to have anomalous dispersion. If Eq. (|74^) is satisfied, the energy conservation 
condition £ l qi +Eq 2 = -E qi + q2 can be satisfied. An excitation Eq 1 can then decay into Eq 1+ q, 2 
by absorbing Eq 2 (Landau), or an excitation Eq 1+ q 2 can decay into two excitations Eq 1 and 
Eq, 2 (Beliaev). 

The condition for such a qi to exist is that the maximum value of |V q .Eq| as a function 
of (q x ,q y ) is greater than c. Due to |V q -E q | q= o = c, we only need to consider the condition 
that |V q -Eq| takes its maximum at q 7^ 0. When |V q -E q | has its maximum value, 

2 Jd 2 sin q x d 

Vq|VqEql = /• 2 =h \ 

'sin q x d + sin q y d 





E a / cosq x d \ 2Jf/ 2 (n c0 ) 2 / . 2 J . 2 JN . , . 

^^(sin'g^ + sin 2 ^^) ( | | = 0. (75) 

cos q y d I -^q 

From Eq. (|75|) . |V q -E q | has its maximum value when cos q x d = cos q y d, i.e., at the two values 
q x = ±q y . If we assume q x = q y and define u = sin 2 q x d, Eq. (|75j) reduces to 



. 2 , -(3a-4) + v / 5a 2 + 24a + 16 

u = sin q x d = — . (76) 

lb 

From Eq. ([76)1 . one sees that u decreases as a increases, vanishing when a = 12. Therefore, 
we conclude that the Landau damping when a < 12 is due to excitations with momentum 
q x = ±q y , and all damping processes in a 2D optical lattice will vanish when a > 12. We have 
confirmed this analytical result by numerically solving the energy conservation condition. 

The condition for the disappearance of Landau damping in a 3D optical lattice can be 
derived by generalizing the procedure described above for a 2D optical lattice. One sees that 
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|V q -Eq| has its maximum when cos q x d = cos q y d = cos q z d, i.e., q x = ±q y = ±g 2 and 

. 2 , -3(q - 2) + V5a 2 + 36a + 36 

sm 2 q x d= ^ '- — . (77) 

One finds that u — > when a — > 18, and damping in a 3D optical lattice vanishes when 
a > 18. As in the 2D case, Landau damping when a < 18 only occurs for an excitation 
with momentum q x = ±q y = ±q z . 

We discuss the energy conservation condition in a 2D optical lattice in detail. We only 
consider the damping of long wavelength phonon q. Using the approximation for the long 
wavelength phonon E n ~ cq and -Eq+ P ~ E^ + V p -E p ■ q, the energy conservation condition 
-E q + -Eq = -Eq+p can be written as 

a E 2 

2 {ll + ) = ^| {^ 2 {Pxd)q 2 x + 2 sin(p x d) sin(p y d)q x q y + sin 2 (p y d)q 2 y ) . (78) 

When g z > and q y — 0, Eq. (|78|) can be solved easily. Defining X = sin 2 2|- and 
Y = sin 2 ^, the solution of Eq. fZHJ) is 



y --( x + l) + lV 8 ^-tx + a < 79 ' 

One can confirm that Eq. ()79j) reduces to the ID result in Eq. (|63|) when Y = 0. 

Equation (|79*jl is plotted in Fig. [7|for several values of a. We see that as a — > 6 the line 
in the {p x ,p y ) plane which satisfies the energy conservation condition shrinks and vanishes 
when a > 6. Therefore, the Landau damping of an excitation with q y = disappears when 
a > 6. po in Fig. 2.8 is given by Eq. (JBHJ). 

For a long wavelength phonon q with q x = q y > 0, we solve the energy conservation 
condition numerically. The solution is shown in Fig. |H1 There is no solution when a > 12 as 
expected. Figures [7| and |H1 clearly show that the threshold value of a for the disappearance 
of damping strongly depends on the direction of q due to the anisotropy of 2D square lattice. 
This result also holds for 3D square lattice. 

As for the ID case, becomes larger than E^ around the threshold value of a in 2D and 
3D optical lattices. In this case, we cannot use the simple Golden Rule expression of the 
Landau damping given by Eq. (JBTj) . We have to extend it using higher order perturbation 
theory in order to calculate the Landau damping when a is close to the threshold value 3] • 
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p d/K 
' x 

FIG. 7: The solution of the energy conservation condition E^ + Ep = E^ +p for a 2D optical lattice, 
with q x > and q y = 0. po is defined in Eq. (|65|) and we note po — > as a — > 6. 



0.8 




FIG. 8: The solution of the energy conservation condition £" q + E p = E Q+P in a 2D optical lattice, 
for q x = q y = O.lir/d. 
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D. Beliaev damping in a ID optical lattice 



We briefly discuss the Beliaev damping of the Bloch-Bogoliubov excitation in this subsec- 
tion. The Beliaev damping is due to spontaneous decay of an excitation into two excitations, 
and thus we need to satify the energy conservation condition £" q = -Eq-p + E p . We focus on 
a ID case in the following. 

In Fig. El the solution of the energy conservation condition for ID optical lattice E q = 



E q _p + E p is shown in a (q,p) plane. As predicted in Sec. IV B[ one finds that the curve of the 
solution shrinks as a increases and vanishes when a > 6, which indicates the disappearance 
of the Beliaev damping for a > 6. 

For a fixed value of a < 6, Beliaev damping is only possible when q is between the 
threshold momenta go and q c shown in Fig. Q. At the threshold momenta go and g c , two 
excitations E q and E q _ p created by the decay of an excitation E q have the same velocity. 
This effect was first predicted by Pitaevskii in 1959 for the phonon excitation in superfluid 
[45]. The two created excitations have the same quasi- momentum g c /2 and energy E qc /2 
at q = q c . At q = q , one of the generated excitations is a phonon having the sound velocity 
c. Therefore, the other one also has the velocity equal to the sound velocity c. 

In addition to Landau damping and Beliaev damping, one also has intercollisional damp- 
ing arising from two body collisions which transfer atoms between the condensate and ther- 



i 



mal cloud at finite temperatures |33l. l34j . Such processes also involve the energy conservation 
condition for three-excitation processes. Thus, the intercollisional damping also disappears 
when a > QD in a D-dimensional optical lattice. 

VI. CONCLUSION 

In conclusion, we have given a detailed treatment of the damping of Bogoliubov excita- 
tions associated with Bose condensates in an optical lattice at finite temperature using the 
tight-binding Bose- Hubbard model, extending the work we reported in Ref. 21 1. 



We have used the Popov approximation in the Bose-Hubbard model to extend the usual 
theory to finite temperatures. As a by-product, we have calculated the number of condensate 
atoms per lattice site as a function of both the temperature and the lattice well depth s. 
These results may be of general interest. 
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FIG. 9: The solution of the energy conservation condition E q = E p + E q _ p for Beliaev damping in 
a ID optical lattice. 

We calculated the damping of Bloch-Bogoliubov excitations in ID, 2D and 3D optical 
lattices. All previous work p, Isj] on this problem only considered dynamical instabilities. 
These studies did not include the dissipation processes we have considered. We find that the 
Bogoliubov-Popov excitation spectrum E n must exhibit "anomalous dispersion" for damping 
processes to occur. This is analogous to the case of phonon damping in superfluid 4 He. In 
the absense of this "bending-up" of the low q spectrum, energy conservation cannot be 
satisfied. As a consequence, excitation damping is absent when a = Un c0 /J > 6D where D 
is the dimension of the optical lattice. 

The first studies [3] of damping of excitations were limited to ID optical lattices along 
the axis of a cigar-shape magnetic trap, but one would need a much tighter magnetic trap 
(in the radial direction) for our ID model results to apply. In the more recent experiments 



by Stoferle et. al. 17[, a 3D optical lattice is prepared first, and the lattice potential depths 
in two lattice axes are then made much larger than the third one to produce 2D array of 
tightly bound ID optical lattices. An analogous 2D optical lattice can be formed by choosing 
a much larger lattice potential depth along one lattice axis than that in other two axes. This 
effectively 2D optical lattice might be better for checking our theoretical predictions than 
the 2D poriodic array of ID tubes used in Ref. [241 , since excitations along the direction 
perpendicular to the 2D lattice potential can be neglected. 
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Due to our use of a tight-binding approximation, our results are not directly applicable 
to the damping of excitations found in very weak optical lattices (s < 1) 18 . Extension of 
our calculations to such weak optical lattices would be clearly of interest |46|. 
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